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Assuming that the hard 7-ray emission of Cen A is a result of synchrotron radiation of ultra- 
relativistic electrons, we derive a lower bound on the local extragalactic magnetic field, B > 10~* 
G. This result is consistent with (and close to) upper bounds on magnetic fields derived from 
consideration of cosmic microwave background distortions and Faraday rotation measurements. 



Increasing our knowledge of the properties of extra- 
galactic magnetic fields (EGM) on scales between Mpc 
and the visible horizon could significantly impact our un- 
derstanding of phenomena in cosmic ray physics (prop- 
agation, dissipation via synchrotron radiation), cosmol- 
ogy (distortion of cosmic microwave background (CMB)), 
and particle physics (early universe primordial magneto- 
genesis during phase transitions). At present, surpris- 
ingly little is actually known about EGM: there are some 
measurements of diffuse radio emission from the bridge 
area between the Coma and Abell superclusters 
which, under assumptions of equipartition allows an es- 
timate of 0(0.2 — 0.6) /iG for the magnetic field in this 
region. Such a large field may possibly be understood if 
the bridge region lies along a filament or sheet ||]. Fara- 
day rotation measurements |^,|^ have thus far served to 
set upper bounds of 0(10^^ — 10^^) G on EGM on var- 
ious scales [DJ^, as have the limits on distortion of the 
CMB The Faraday rotation measurements sample 

EGM of any origin out to quasar distances, while the 
CMB analyses set limits on primordial magnetic fields 
l^. We will discuss these bounds in some detail in the 
concluding section, when comparing with our results. 

In this paper, we will use the results of a set of cosmic 
ray measurements taken over a ten-year period at the 
Sydney University Giant Air-shower Recorder (SUGAR) 
pO| , and covering much of the southern sky, in order 
to set a lower bound on the size of the i?-field between 
Earth and the nearby galaxy Centaurus A (Cen A), at 
a distance d ~ 3.4 Mpc. The bound will be based on 
current models of dynamics in the "hot spots" of radio- 
galaxies, and these assumptions will be tested by future 
observations of ultra-high energy cosmic ray (CR) and 
high energy gamma ray fluxes from the direction of Cen 
A. We will initially outline the dynamical picture con- 
cerning the generation of high energy gamma rays. This 
picture leads to our Eq. (|^), and then as a consequence 
to our lower bound on B. We may immediately note that 
we are aware of the crude angular and energy resolutions 
of the SUGAR observations. Our analysis will take into 
account these uncertainties. 

Fanaroff-Riley II (FRII) galaxies are the largest 
known dissipative objects (non-thermal sources) in the 
universe. Localized regions of intense synchrotron emis- 
sion, so-called "hot spots", are observed within their 
lobes. These regions are presumably produced when the 



bulk kinetic energy of the jets ejected by a central active 
nucleus (supermassive black hole -I- accretion disk) is re- 
converted into relativistic particles and turbulent fields 
at a "working surface" in the head of the jets In 
what follows, we will adopt the first order Fermi shock 
acceleration mechanism in hot spots of FRII glaxies, as 
discussed in Refs. [|l3|,0, to account for particle accel- 
eration to ultra high energies The apphcability of 
this scenario to Cen A will be discussed in the next para- 
graph. The subtleties surrounding the conversion of par- 
ticles' kinetic energy into radiation provide ample mate- 
rial for discussion [^6|-Q . The most popular mechanism 
to date relates 7-ray emission to the development of elec- 
tromagnetic cascades triggered by secondary photomeson 
products that cool instanstaneously via synchrotron ra- 
diation p6| , p7t . The characteristic single photon energy 
in synchrotron radiation emitted by an electron is 

For a proton this number is {mp/rrif,)^ ~ 6 x 10^ times 
smaller. Here, B^g is the magnetic field in units of /iG 
and Eig = eV. Thus, it is evident that high en- 

ergy gamma ray production through proton synchrotron 
radiation requires very large (0(100 G)) magnetic fields. 

We now discuss how the less luminous FRI galaxy Cen 
A fits into this picture. This radio-loud source {I « 310°, 
b ~ 20°), identified at optical frequencies with the galaxy 
NGC 5128, is the closest example of the class of active 
galaxies . Different multi- wavelength studies have re- 
vealed a rather complex morphology: it comprises a com- 
pact core, a jet also visible at AT-ray frequencies, a weak 
counterjet, two inner lobes, a kpc-scale middle lobe, and 
two giant outer lobes. The jet would be responsible for 
the formation of the northern inner and middle lobes 
when interacting with the interstellar and intergalactic 
medium, respectively. There appears to be a compact 
structure in northern lobe, at the extrapolated end of 
the jet. This structure resembles the hot spots like those 
existing at the extremities of FRII galaxies. However, at 
Cen A it lies at the side of the lobe rather than at the 
most distant northern edge, and the brightness contrast 
(hot spot to lobe) is not as extreme [ pO[ . 

In order to ascertain the capability of Cen A to ac- 
celerate particles to ultra high energies, one first applies 
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the Hillas criterion |Q for locaUzing the Fermi engine 
in space, namely that the gyroradius Tg = IIOE20/ B^q 
kpc (£'20 = -E/IO^" eV) be less than the size of the mag- 
netic region. Low resolution polarization measurements 
in the region of the suspected hot spot give fields as high 
as 25 /iG j20|. In certain of the regions where measure- 
ments at both high and low resolution are available, the 
_B-field at high resolution can be twice that at low reso- 
lution. The higher resolution can reveal amplification in 
the post-shock regio n [^21 , yielding B- fields possibly as 
high as 50 — 60 /iG |20y23| ] . The radio- visible size of the 
hot spot can be directly measured from the large scale 
map of Ref. Q, giving i?HS — 2 kpc. The actual size 
can be larger by a factor ^ 2 because of uncertainties 
in the angular projection of this region along the line of 
sight |2^. If one assumes that the magnetic field of the 
hot spot is confined to the visible region, then the limit- 
ing energy is ~ 2 x 10^" eV. However, it is plausible that 
the shock structure in the hot spot extends beyond the 
radio- visible region | p^ , p6[ . 

In light of this, we apply to Gen A the analyses of 
Refs. in which a limiting energy is obtained by 

balancing the characteristic time scale for diffusive shock 
acceleration 



kpc. Since the range of values for a and the jet velocity 
conform to expected values, it is plausible that Gen A 
can accelerate particles to energies > 10^° eV. 
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FIG. 1. Jet velocity as a function of parameter a (defined 
in text), for different proton energies. 
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In the above, pjet = jet velocity in units of c, u is the ratio 
of turbulent to ambient magnetic energy density in the 
region of the shock (of radius R), B is the total magnetic 
field strength, a is the ratio of photon to magnetic energy 
densities, ctt is the classical Thomson cross section, and 
A is a measure of the relative strength of jp interactions 
against the synchrotron emission. In Ref. [Q, A is es- 
timated to be « 200 almost independent of source. In 
Eq.(H) we assume that the turbulent component of the 
magnetic field follows a Kolmogorov spectrum with spec- 
tral index 5/3. Equating these characteristic times yields 
a value for the maximum proton energy 



E20 = 1.4 X 10^ 
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where i?kpc = kpc. One can estimate u 0.4 from 
the radio spectral index of synchrotron emission in the 
hot spot and the observed degree of linear polarization in 
the same region . The jet velocity is model dependent, 
and estimates range from ~ 500 km s~^ to 0.99c [Q. For 
FRI galaxies, a is expected to be ^ 1 ||l^], and in our 
analysis we will sample a region of small a. In Fig. 1 we 
plot the relation between /3jet and a required to attain 
various energies, for fiducial values B = 60 /xG, i? = 4 



Recent observations of the gamma ray flux for ener- 
gies > 100 MeV (by the Energetic Gamma Ray Ex- 
periment Telescope (EGRET) [M) allow an estimate 



10 erg s for the source 



This value of 



is consistent with an earlier observation of photons in 
the TeV-range during a period of elevated activity 
and is considerably smaller than the estimated bolomet- 
ric luminosity iboi ~ lO^'^erg s"^ ||l^. Data across the 
entire gamma ray bandwidth of Gen A is given in Ref. 

reaching energies as high as 150 TeV Data 
at this energy await confirmation. For values of B in 
the fiG range, substantial proton synchrotron cooling is 
suppressed, allowing production of high energy electrons 
through photomcson processes. The average energy of 



synchrotron photons scales as E^ ~ 0.29Ej 



With 



this in mind, it is straightforward to see that to account 
for TeV photons Gen A should harbor a population of 
ultra-relativistic electrons with E 6 x 10^^ eV. We fur- 
ther note that this would require the presence of protons 
with energies between one and two orders of magnitude 
larger, since the electrons are produced as secondaries 

There are plausible physical arguments ||17|,^ as well 
as some observational reasons to believe that when 
proton acceleration is being limited by energy losses, the 
GR luminosity Lck ~ L^. In the spirit of |37|, we intro- 
duce e, the efficiency of ultra high energy GR production 
compared to high energy 7 production — from the above, 
we expect e ~ 1. Using equal power per decade over the 
interval (-Emin, ^-max), we estimate a source luminosity 
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where L41 = luminosity of Cen A/lO^^erg s"-'^ and the 
subscript "0" refers to quantities at the source. 

For a relatively close source like Cen A, one can ne- 
glect interactions of cosmic ray protons with the univer- 
sal background radiations. Consequently, the shape of 
the spectrum would be unmodified . In order to dis- 
cuss directionality, we first consider the case where there 
is no intergalactic magnetic field. Then en route to us 
the protons suffer no deflections and no spectral distor- 
tions. With the source luminosity (||), the counting rate 
at Earth in an energy bin E^) for a detector of area 
S due to Cen A would be 
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(Here Ei^20 = Ei/lO'^^eV, etc.) This would all be con- 
centrated in a cone of half- angle ^res, the observational 
resolution, about the direction of Cen A. 

As mentioned previously, the only existing measure- 
ments of ultrahigh energy cosmic rays arriving from di- 
rections which include Cen A were carried out at SUGAR 
during a ten year period 1968-1979. The array was lo- 
cated at 30°32' S, 149°46' E, and presented a total area 
of 100 km^. Shown as asterisks in Fig. 2, as given in Ref. 
pof , are the arrival directions of 80 events with energies 
above 4 x 10^^ eV. The direction of Cen A is indicated 
with a five-pointed star, and the experimental la uncer- 
tainties in arrival directions of events near the Cen A 
direction are indicated by dashed lines |3^. The solid 
line ovals surrounding Cen A indicate regions within 10° 
and 25°, respectively, of Cen A, and their significance 
will be explained shortly. 

We now obtain the B = event rate expected at 
SUGAR in the direction of Cen A, as predicted by 
Eq.(^. The energy bin is appropriate to the data of 
pof , -El, 20 = 0.4, E2.2Q — 2. A conservative lower bound 
on the energy of protons which are progenitors of high 
energy gamma radiation is £^min = 1 x 10^^ eV. It also 
reasonable to take for the upper bound Emax = 4 x 10^° 
eV Q. With this input, we obtain from Eq.(||) 
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In the preceding discussion we have indicated that under 
plausible conditions may expect e ^ 1. However, even if 
eL4i is as small as 0.1, we may, in the absence of a mag- 
netic field, expect 8 events in 10 years from the direction 
of Cen A. 

The cosmic ray orbits undergo bending in both the 
Galactic and extragalactic magnetic fields. The magnetic 
deflection of protons in the Galactic disk has been stud- 
ied in detail in [O. This analysis includes two extreme 



options for the behavior of the field, reflecting the differ- 
ent symmetries with respect to field reversals in the r- 
and z-directions. The B-fleld has a l/r behavior, with 
deviations calculated out to 20 kpc from the Galactic 
center. The r.m.s. deviation (averaged over arrival di- 
rection), for an energy of 4 x 10^^ eV, varies from 7.9° to 
10.5° in going between the two models. This deviation 
shows an approximately linear decrease with increasing 
energy. All the events which we will consider have ener- 
gies > 4 X 10^^ eV, so that conservatively, we allow an 
uncertainty of 10° in arrival direction due to deviation 
in the Galactic magnetic field. Events within the inner 
oval of Fig. 2 represent events which could have origi- 
nated from Cen A, and have suffered deflection only in 
the Galactic field. It is clear that at most 1 event can be 
categorized in this manner, so that in order to have at 
least 8 events (see previous paragraph), some amount of 
extragalactic field is necessary. 

We treat the extragalactic deviation in a standard 
manner 0. With B^g = -B/IO^^G, the Larmor ra- 
dius of a particle in this field is tl — 10^ Mpc i?2o/SnG- 
If this is sizeably larger than the coherence length of the 
magnetic field €coh, the typical defiection angle from the 
direction of the source, located at a distance d, can be es- 
timated assuming that the particle makes a random walk 
in the magnetic field [M 
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Thus, in the case that 0{E) <^ 1, we expect that the 
counting rate (|^) would be spread out over a cone of 
half-angle 



9{E) ~ {e{Ef 
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We now draw a second oval in Fig. 2 to include the 
least number of events compatible with the expectation 
from Cen A. This second oval contains (or partially con- 
tains) 7 events within 25° of Cen A. From our previous 
discussion of the Galactic effects, this leads us to con- 
clude that cosmic rays with E > A x 10^^ eV experience 
a deviation of at least 15° in extragalactic magnetic fields 
during their transit from Cen A to Earth. 

In conjunction with Eqs.(^ and (^ we utilize this con- 
clusion to obtain a lower bound on the _B-field between 
Earth and Cen A. The average energy E of the cosmic ray 
events above 4 x 10^^ eV, as reported in |l^, is 6.3 x 10^^ 
eV. Requiring e(E) > 15°, with 6'rcs = 3° and d = 3.4 
Mpc, we obtain from (ph and (0) 



B>9.5x IQ-^ (4oh/lMpc) 
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Since the coherent length is most likely < 1 Mpc, we 
obtain the bound stated in the abstract, 



B >1.0x IQ-^ G 



(11) 
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FIG. 2. The nominal arrival directions (a = right ascen- 
sion, S — declination) of SUGAR events with energies above 
4 X 10^^ eV. Also shown in solid lines are contour maps in- 
dicating the circular areas of the celestial sphere centered at 
Cen A (indicated by ★) with 10° and 25° radii. The dashed 
lines surrounding several of the events indicate the angular 
resolution of the experiment. 

Our analysis has utilized the SUGAR data with en- 
ergies assigned according to the Hillas model |43). With 
this prescription, the energy spectrum agrees with that of 
other experiments An alternate method [Q which 
underestimates the energy by a factor of more than 2 
leads to a spectrum which does not show this agreement 
pOf . The lower energy would have two effects in our anal- 
ysis: [i) the expected event rate (0) would double because 
of the smaller lower limit on the energy bin (m) the devia- 
tion in the Galactic magnetic field would roughly double 
leading to an increase in the size of the inner oval 
in Fig. 2. Following the same procedure as previously, 
and still utilizing the conservative bound eL^i > 0.1, 
we would find that B > 5 nG. However, because there is 
general consensus |^ that the spectral evidence strongly 
favors the Hillas energy normalization, we maintain our 
bound as S > 10 nG le). 

The limit ( pT| ) has an immediate implication with re- 
spect to lower energy cosmic rays: from Eq.(||), one can 
conclude that ii B > 10 nG all directionality is lost for 
protons with energies below ~ 2 x 10^^ cV. This is con- 
sistent with the absence of any anisotropy in the events 
observed at SUGAR at these lower energies ]47| . 

We turn to discuss our result in light of considera- 
tions which impose upper bounds on magnetic fields. (1) 
An analysis of He^ production provides a weak bound 
Bq < 10~^ G on a homogeneous primordial field p8[ ; 
however, such fields are largely dissipated prior to nu- 
cleosynthesis (2) From limits on distortions of the 
cosmic microwave background (CMB), there is an up- 



per bound Bq < 3.4 x 10~ G for fields homogeneous on 
present horizon scales . Our lower bound (^l]) exceeds 
this, so that we may conclude that the magnetic field 
from here to Cen A contains components in excess of a 
homogeneous relic background field. (3) The analysis of 
CMB distortions extended to inhomogeneous fields pro- 
vides a larger bound Bq < 3 x 10~* G on scales between 
400 pc and 0.6 Mpc This hes above our limit, so 
that the local EGM is consistent with being of primor- 
dial origin. (4) The absence of a positive signal in Fara- 
day rotation measurements on QSO's |q| provides upper 
limits on magnetic fields (of any origin) as a function of 
reversal scale. These bounds depend significantly on as- 
sumptions about the electron density profile as a function 
of red shift z. When electron densities follow that of the 
Lyman-a forest, the average magnitude of the magnetic 
field receives an upper limit oi B < 10^^ G for reversals 
on the scale of the horizon, and B < 10^* G for reversal 
scales on the order of 1 Mpc |^ . The latter upper bound 
is roughly coincident with our lower limit in our galactic 
neighborhood. Local fluctuations in electron densities on 
scales of 1-30 Mpc can lead to very large concommitant 
fluctuations in the corresponding magnetic fleld so 
that the bound in ||] should be read as averaged over 
many reversal cells between Earth and the light sources 
at distances out to z = 2.5. If it should happen that the 
fluctuations in the magnitude of B in the Mpc cells out to 
the horizon are of the same order as i?avg itself, then our 
result can imply that the average value of the intergalac- 
tic magnetic field is of the order of 10~® G. Because of 
the B^ dependence of synchrotron radiation loss by sec- 
ondary electrons, a fleld of this strength can have impor- 
tant implication on the development of electromagnetic 
cascades. Such cascades are characteristically associated 
with the decay of supermassive relic particles or topolog- 
ical defects. An average magnetic field of 10~^ G will 
impose strong constraints on "top-down" models |]50|] . 

The SUGAR observations used in this analysis were 
recorded more than twenty years ago. Both resolution 
and statistics will be vastly improved with data to be 
available from facilities which will observe the southern 
sky. These are now coming on line (Auger [pl|), or are 
expected to (EUSO / OWL ) in the not-too-distant fu- 
ture. These data will significantly enhance our knowledge 
of magnetic field strengths in the extragalactic neighbor- 
hood of the Milky Way. 
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